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Abstract 

H 

C^ A homogeneous Poisson-Voronoi tessellation of intensity 7 is observed in a convex body W. We associate to 

^^ each cell of the tessellation two characteristic radii: the inradius, i.e. the radius of the largest ball centered at 

^^ the nucleus and included in the cell, and the circumscribed radius, i.e. the radius of the smallest ball centered 

T-H at the nucleus and containing the cell. We investigate the maximum and minimum of these two radii over all 

C^ cells with nucleus in W. We prove that when 7 — >■ 00, these four quantities converge to Gumbel or WeibuU 

distributions up to a rescaling. Moreover, the contribution of boundary cells is shown to be negligible. Such 

approach is motivated by the analysis of the global regularity of the tessellation. In particular, consequences of 
our study include the convergence to the simplex shape of the cell with smallest circumscribed radius and an 
upper-bound for the Hausdorff distance between W and its so-called Poisson-Voronoi approximation. 
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1 Introduction 

^D Let X be a locally finite subset of R'* endowed with its natural norm | • | . The Voronoi cell of nucleus a; G x is the 

■^ set 

cn C^i^) -{ye R^ \y-x\<\y- x'lx ^x' e x}- 

ll" When X — ^j is a homogeneous Poisson point process of intensity 7, the family {Cx^{x),x e X^} is the so- 

. ^H called Poisson-Voronoi tessellation. Such model is extensively used in many domains such as cellular biology [32] , 

JC astrophysics [22|i telecommunications [51 and ecology [^. For a complete account, we refer to the books i3D], [37], 

H [27] and the survey [B]. 

To describe the mean behaviour of the tessellation, the notion of typical cell is introduced. The distribution of 
this random polytope can be defined as 



^[^(''^-''i - ^,^ 



E /(Cx,(i)-i) 

xex^nB 
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where / : J(^'^ — > IR is any bounded measurable function on the set of convex bodies Jif^ (endowed with the 
HausdorfF topology), A^ is the d-dimensional Lebesgue measure and B is a Borel subset of R'* with finite volume 
Xd{B) £ (0,00). Equivalently, "^^ is the Voronoi cell Cx u{o}(0) when we add the origin to the Poisson point 
process: this fact is a consequence of Slivnyak's Theorem, see e.g. Theorem 3.3.5 in [37]. The study of the typical 
cell in the literature includes mean values calculations J25] , second order properties [TT] and distributional estimates 
[5], |3], pS]. A long standing conjecture due to D.G. Kendall about the asymptotic shape of large typical cell is 
proved in [15] . 

To the best of our knowledge, extremes of geometric characteristics of the cells, as opposed to their means, have 
not been studied in the literature up to now. In this paper, we are interested in the following problem: only a part 
of the tessellation is observed in a convex body W (i.e. a convex compact set with non-empty interior) of volume 
Xd{W) — 1 where A^ denotes the Lebesgue measure in R'*. Let / : J^^^ — > R be a measurable function, e.g. the 
volume or the diameter of the cells. What is the limit behaviour of 

M/(7)= max fiCx^ix)) 

when 7 goes to infinity? By scaling invariance of X.y , it is the same as considering a tessellation with fixed intensity 
and observed in a window Wp := pW with p — >■ cxd. We give below some applications of such approach. 

First, the study of extremes describes the regularity of the tessellation. For instance, in finite element method, 
the quality of the approximation depends on some consistency measurements over the partition, see e.g. [TH] . 

Another potential application field is statistics of point processes. The key idea would be to identify a point 
process from the extremes of its underlying Voronoi tessellation. A lot of inference methods have been developed 
for spatial point processes [IH]. A comparison based on Voronoi extremes may or may not provide stronger results. 
At least, the regularity seems to discriminate to some extent to some point processes (see for instance a comparison 
between a determinantal point process and a Poisson point process in |23jV 

A third application is the so-called Poisson- Voronoi approximation i.e. a discretization of a convex body W by 
the following union of Voronoi cells 

yx,{w)^ U cxM- 

xexnw 

The first breakthrough is due to Heveling and Reitzner [14j and includes variance estimates of the volume of 
symmetric difference. However, the Hausdorff distance between the convex body and its approximation has not 
been studied yet. It is strongly connected to the maximum of the diameter of the cells which intersect the boundary 
of dW . We discuss this in section |4] and prove a rate of convergence of the approximation to the convex body with 
a suitable assumption on W . 

Concretely, we are looking for two parameters af{'j) and ^7(7) such that 

a/(7)M/(7) + ^^/(7) A Y 

7—^00 

where K is a non degenerate random variable and — > denotes the convergence in distribution. Up to a normaliza- 
tion, the extreme distributions of real random variables which are iid or with a mixing property are of three types: 
Frechet, Gumbel or WeibuU (see e.g. [21] and [2T|). More about extreme value theory can be found in the reference 
books by De Haan & Ferreira [5] and by Resnick [33]. Some extremes have been studied in stochastic geometry, 
for instance the maximum and minimum of inter-point distances of some point processes (see |31j . [25] and \W\ ) 
or the extremes of particular random fields [2D] but, to the best of our knowledge, nothing has been done for 
random tessellations. In our framework, the general theory cannot directly be applied for several reasons: unknown 
distribution of the characteristic for one fixed cell, dependency between cells and boundary effects. Moreover, the 
exceedances can be realized in clusters. For example, when the distance between the boundary of the cell and its 
nucleus is small, this is the same for one of its neighbors. Such clusters lead to the notion of extremal index, which 
was introduced by Leadbetter in [22J, and that we will study in a future work. 



In this paper, we are interested in the characteristic radii i.e. inscribed and circumscribed radii of the Voronoi 
cell Cx {x) defined as 

r{Cx,{x)) = max{r > 0, B(x, r) C Cx,{x)} and R{Cx,{x)) = min{i? > 0,B(x,i?) D Cx,{x)} 

where B{x, r) is the ball of radius r centered at x. Two reasons led us to the study of these quantities. First, the 
distribution tails of the inradius and circumscribed radius of the typical cell are easier to deal with [4J compared to 
other characteristics such as the volume or the number of hyperfaces. Secondly, knowing these two radii provides a 
better understanding of the cell shape since the boundary of Cx {x) is included in the annulus B{x, R{Cx {x))) — 
B{x,r{Cx (x)))- We consider the extremes 

rmax(7) = max ''(Cx^(x)), rmi„(7) = "^im r{Cx-,{x)) 

xex.,nw xeXyHW 

i?max(7) = max i?(Cx.,(a:)), i?min(7) = mm i?(Cx.,(a:)). 

xEX^nw xex^nw 

In the following theorem, we derive the convergence in distribution of these quantities over cells with nucleus in W . 
Theorem 1. Let X^ be a Poisson point process of intensity 7 and W a convex body of volume 1 in R'*. Then 

P {2''Kdir^,Ulf - log(7) <t) ^ e~''\ t e R, (2a) 

7— J'OO 



{2''-^Kd7^rn,inhf >t) ^ e-\ t > 0, (2b) 

7— >oo 



P {iidlRn...hf - log (ai7(log7)'*-i) <t) ^ e~'~\ t e R, (2c) 



P(a2 W'^+')/('^+^)i?.„in(7)' > i) ^ e-*'"\ t>0, (2d) 

7—^00 



where ai and a2 are given in (43 1 and (17 1 and Hd = '^^(^(0, 1)) 



The limit distributions are of type II and III and do not depend on the shape of W. One can note that the ratios 
^max(7)/''min(7) and i?max(7)/-Rmin(7) are of rcspcctivc orders (7log7)^/'' and [^^^('^+^^ log7)^/'^. This quantifies to 
some extent the irregularity of the Poisson- Voronoi tessellation. Moreover, the ratio ?'i„ax(7)/^max(7) is bounded. 
It suggests that large cells tend to be spherical around the nucleus. This fact seems to confirm the D.G. Kendall's 
conjecture. 

As it is written. Theorem 1 is not applicable for concrete data. Indeed, in practice, the only cells which can be 
measured are included in the window. The following proposition addresses this problem. 

Proposition 2. The extremes of characteristic radii over all cells included in W or over all cells intersecting dW 
have the same limit distributions as r-a_-iax{l), ?'min(7), ^max(7) o^nd Rmin{l)- 

The convergences are illustrated in Figure^for the cells which are included inW = [0, 1]^. For sake of simplicity, 
the Poisson point process has been realized only in W . Because of Proposition [2] and related arguments, this does 
not affect the distribution over cells included in W. Simulations suggest that the rates of convergence are not the 
same for all these quantities. Indeed, in a future work, we will show that the rate is of the order of 7~^, ^~^/'^ and 
^-1/6 for rmin(7)> ''max (7) and Rminh) respectively. 
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Figure 1: Empirical densities of the extremes based on 3500 simulations of PVT in 2D with 7 = 10000, for the cells 
included iaW— [0, 1]^, on Matlab®. (a) Cell maximizing the inradius. (h) Cell minimizing the inradius. (c) Cell 
maximizing the circumradius. (d) Cell minimizing the circumradius. 



All results of Theorem [T] use geometric interpretations. For the circumscribed radii i?max(7) and i?min(7), we 
write the distributions as covering probabilities of spheres. The inscribed radii can be interpreted as interpoint 
distances. A study of the extremes of these distances has been done in several works such as [W and [13]. For 
sake of completeness, we have rewritten these results in our setting in particular because the boundary effects are 



highly non trivial. Convergences ( |2a[ ) and ( 2d ) could be obtained by considering underlying random fields and using 
methods inherited for [Tj and [Ml. However, this approach does not provide (2b I and (2c I. We will develop this 



idea in a future work and deduce some rates of convergence. 

The paper is organized as follows. In section 2, we provide some preliminary result which shows that the 
boundary cells are negligible and implies Proposition [2] In sections 3, 4 and 5, proofs of (2d I, (2a I, (2c I and (2b I 
are respectively given. Section 3 requires a technical lemma about deterministic covering of the sphere by caps 



which is proved in appendix. Section 4 contains an application of ( 2c I to the Hausdorff distance between W and its 



Poisson-Voronoi approximation. In section 5, we get a specific treatment of boundary effects which is more precise 
than in section 2. 

In the rest of the paper, c denotes a generic constant which does not depend on 7 but may depend on other 
quantities. The term u^ denotes a generic function of t, dependending on 7, which is specified at the beginning of 
sections 3,4 and 5. 



2 Preliminaries on boundary effects 

In this section, we show that the asymptotic behaviour of an extreme is in general not affected by boundary cells. 
We apply that result directly to the extremes of characteristic radii in order to show that Theorem [T] implies 
Proposition [2J 

Let / : J(fd H> IR be a fc-homogeneous measurable function, Q < k < d (i.e. f{XC) — X^ J{C) for all A G 1R+ and 
C £ Jtd)- We consider for any I G R 



M'AjJ) 



max f{Cx^{x)), 



Mfij,l)^ max nCx,{x)), 

M}(7,0= , max f{Cx^{x)), 

where Wi+i = (1 + l)W. When I = 0, these maxima are simply denoted by M5(7), Mf{'y) and MU'f). We define, 
for all e > 0, a function Ij as 

l^ = j-(^-^)/d, (4) 

Under suitable conditions, the following proposition shows that Af5(7), Mf{'y) and MUj) satisfy the same conver- 
gence in distribution. 

Proposition 3. Let Y be a random variable and a^, b^ two functions such that 

""^ ^ 1, l^b^ -^ and ^"7± - «7^± _^ Q (5-) 



a-y^ 7->-oo 7->-oo a-y 7->oo 

with 7-|_ = (1 + ^)'^7 and 7_ = (1 — ^)'^7 /or a certain e. Th 



en 



a^M'l{j)+b^ A Y ■i=^a^Mf{j)+b^ A y -^^^ a^MU7) + 6^ A K 

^ 7^oo 7— >oo ■' 7— >oo 

Before proving Proposition^ we need an intermediary result due to Heinrich, Schmidt and Schmidt (Lemma 4.1 
of [12]) which shows that, with high probability, the cells which intersect dW have nucleus close to dW. Actually, 
they showed it for any stationary tessellation of intensity 1 which is observed in a window pW with p — > 00. For 
sake of completeness, we rewrite their result in a more explicit version for a Poisson-Voronoi tessellation. 

Lemma 1. (Heinrich, Schmidt and Schmidt) Let us denote by A^ and B^ the events 

A = < n {Cx,{x) nW^0}U{xe Wi+i^} I and B^ = J f| {Cx^{x) C W} U {x <^ W^i_/J 

where l^ is given in Q. Then P{Aj) and P{By) converge to 1 as 7 goes to infinity. 
Proof of Lemma 111 In [12] , it is shown that 



fl {Cx, (x) nWp = 0}u {x e W^p+,(p)} [ n I n {Cx, {x) c W^} u{x^ W^p-g(p)} I 



^ 1 (6) 

7— >-oo 



where q{p) is the solution of the functional equation 

/ = i7(<Z^(p)). 

The function H : IR+ — > IR_|_ is convex, strictly increasing on its support {xo, 00) (for some xq > 0), such that H{x) /x 
is non-decreasing for x > 0, \miH{x)/x = 00 and E[i/(£''^('ifi))] < cx3 where -D('^i) is the diameter of the typical 
cell. 

In the case of a Poisson-Voronoi tessellation, q{p) can be made explicit. Indeed, we can show that all moments 
of D{'^i) exist since I?(^i) < 2R{^i) and -R('^i) is shown to have an exponentially decreasing tail in any dimension 
by an argument similar to Lemma 1 of [3]. Consequently, for a fixed e € (0, 1), the functions H and q can be chosen 
as H{x) = x^l'^ and (/(p) = p"^ . Using the scaling property of Poisson point process, 

Xi n w^p i 7'^''(^7 n w) and Xi n w^^^^^p^ = -i^^\x^ n i^i±^) 

5 



where ^ ~ p'^ and l^ ~ ^ '-^ '^•'/'^. We deduce Lemma 111 from (l6|. D 

Proof of Proposition pi First equivalence: Let us assume that a^M^A'^) + h^ converges in distribution to Y. On 

the event A^,'ix ^ X^, ^x^{x) nW ^ =^ x e VFi+^. Hence 

M){i)<Mf(nA^)<M)(j,l^). (7) 

Because of Lemma [l] it is enough to show the convergence in distribution of the random variables conditionally on 
A^. Thanks to the scaling property of Poisson point process and the /c-homogeneity of / 

Af)(7,M = (l + ^)'^^/(7+) (8) 

with 7+ = (1 + I'y)''^- We deduce from (l7|,(|8]) and (Is]) that a^Mfl-j, l-y) + b~^ converges in distribution to Y . By the 
scaling property, we get 

a^Mf{-f)+h^ A Y. (9) 

7—^00 

Conversely, if dol) holds then, using the fact that 

and proceeding along the same lines, we get a^M5(7) + 6^ — > Y . 

-' 7— >oo 

Second equivalence: On the event B-y, \/x £ X^, x G Wi-i =^ Cx {x) C W . We prove the second equivalence 
as previously noting that, conditionally on B^ 

M}(7, -ly) < Mf{^, -I,) < M}(7) < Mf{j). (10) 

D 



Mf{j) < Af)(7) < Mf{j,l,) I (1 + ly)Mf{^+) 



3 Proof of (2d) and (2a) 



Proof of (2d I. Let t > be fixed. We denote by Uj the following function: 

l/d 



rit) - {a^'n-/j'^'^+'y^^^^h) (11) 



where a2 is given by (17). Our aim is to prove that IP(i?min(7) ^ u-y) converges to e * where i?min(7) has been 
defined in (fTl). The main idea is to deduce the asymptotic behaviour of i?min(7) from the study of finite dimensional 
distributions (-R(Cx-,u{xk}(^i))' • • ■ ; ^(Cx-,u{xk}(2:a:)) for all {x^} = {xi...,xk} and K > 1. To do this, we 
write a new adapted version of a lemma due to Henze (see Lemma p. 345 in [13') in a context of Poisson point 
process. 

Lemma 2. Let / : Jtd — > R, i^ : J^ -> R be two measurable functions and A a Borel subset of IR. Let us assume 
that for any K > 1, 

7^ / P (V^ < KJ{Cx^u{^^}ix,)) < %, i^(Cjf^u{x^}(2:.)) G A) d^K -^ A^ (12) 

where dx^ = dxi ■ ■ ■ dxK- Then 



f{Cx,{x))>u-A ^ e 

x£X^r\W,F(Cx^(x))£A ^ ' j I^OQ 



mm 



-\ 



Proof of Lemma [2} Let K be a fixed integer. The proof is close to the proof of Henze's Lemma and uses 
Bonferroni inequahties: one can show that if A^^x is an X-^-measurable event for all x G X^ n W , then 



2-^ ( lU-+l 



k=0 



k\ 



E 



1a. 



■1a. 



(xi,...,Xk)^<^X^r]W 



<P| U A.,x. 
,xex-,r\w 






fe=0 



fc! 



(xi,...,a;fc)^eX^.nH' 



(13) 



where (xi, . . . , Xfc)^ means that (xi, . . . , x^) is a fc-tuple of distinct points. Applying ( 13 1 to 

Ax,x-, = {I{Cx,{x)) < u^} n {F{Cx,{x)) e A}, 
from Slivnyak's formula (see Corollary 3.2.3 in [37]), we obtain 

E S:i^^' /..., P (^* ^ ^,/(C^Jf-,u{x..}(a:.)) < %,F(Cx-,u{x^}(a;»)) e A) dx^ 



fe=0 

< p 



fc! 



W' 



min /(Cv (x)) > u^ 

xex^nw.F{Cx.Ax))eA ' 



2K 

^E 

fc=0 



^^7^ / P (Vz < K,J{Cx,u{.^}{x^)) < u„F{Cx,u{.^}{x,)) E A) dx,. 

It- Jw'' 



From (12), we obtain 



2K+1 



Y^ ^^X'' < liminf P I mm 



fc=0 



fc! 



7^oo \xex^,nw,F{Cx^{x))eA 



f{CxJx))>u, 



2K 



< lim sup P 

^-yoo \xex.,nw,F{Cx^{x))eA 



min /(Cx,(x))>^,) <5]^ ^^'^'= 



A:=0 



fc! 



-A^ 



We conclude the proof by taking K ^ oo. 



D 



We apply Lemma[2]to f{Cx^{.x)) — R{Cx-,ix)). The function F{Cx^{x)) — Fd^i{Cx^{x)) denotes the number 
of hyperfaces of the cell Cx (x). In all the proof, the event considered is ^ = P. We notice that the choice of the 
function F is of no importance here but will be essential in the proof of Propositions |4] and [5] From Lemma |2] it 
is sufficient to study the limit behaviour of 



7^ / P (Vi < K,R{Cxu{.K}ix^)) < %) d^K 
for all integer K. We divide the proof into two parts. 



(14) 



Step 1 When K — \, using the stationarity of X^ and the fact that Xd{W) = 1, we show that the integral ( [14| 
is ^V{R{Cx u{o}(0)) < u^)- As in [B] section 5.2.3, we can reinterpret the distribution function of R{Cx u{o}(0)) 
as a covering probability to get 



Y-.-2''..W,(2'«'^7^^7)' 



7P(i?(Cx,u{o}(0)) < %) = 7 > e-' -^■'^-' 



k=0 



fc! 



Pk 



(15) 



where pk is the probabiHty to cover the unit sphere with k independent spherical caps such that their norniahzed 
radii are distributed as di'{9) — dn sm{TT9) cos'^~^ {TT9)l[Qi/2]{O)d0. The equahty conies from the fact that 

R{Cx u{o}(0)) < "7 ^=^ the family {.s/y{Q),y e X^} covers S{Q,u^f) 

^=> the family {^^(O), y e X^ n i?(0, 2u^)} covers S{0, Uj) 

where 

£/y{x) = S{x, u^) n H+{x) (16) 

and H^(x) is the half-space which contains y and delimited by the bisecting hyperplane of [x,?/]. 




Figure 2: Interpretation of the circumscribed radius as a covering of sphere. 



We denote by 



a2 



2d{d+l) 

JdTi) 



-Pd+i 



i/id+i) 



>o. 



For example, when d = 2, a2 — (^ — :^) 



Since pk = for all k < d, ( |15| gives 



(2^-.7<)'^+\_,..,,„.^^^^^^ g ^-.^....g (2'^..7<)^ 



k=d+2 



k\ 



Pk- 



(17) 



The first term converges to f^^^ from (11 1 and (17). The second term is negligible since 7(7^^^)'^+^ — c- ^ i/Cc'+i) 
converges to as 7 tends to infinity. This shows that 



7 / P(i?(Cx M{^}(x)) < u^)dx -^ t 

JW 7->-oo 



d+1 



(18) 



step 2 When K >2,we use the same interpretation as in step 1: for all xj^ ~ (.ti, . . . , xk) S W^ , and i < K 

RiCx^u{xK}i^i)) < "7 ^=^ tJ^c family {£/y{xi),y £ X^U {x^} - {xi}} covers S^x^.u^) 

<=> the family {^^(a;^),?; G (X-^ U {x/^-} — {xi}) B{xi,2u^^)} covers S{xi,u^). 

Hence, writing the previous event as ^^S{xi,Uj) covered", we have 

P (Vi < K,R{Cx.,u{.^}{x,)) <u^)^pif] {Six,,u^) covered} . (19) 

We have now to consider the spherical caps induced by both the points Xj,j 7^ i and the points from Xj. For 
all Xk — {xi, . . . , Xk) G W^ , we denote by nj(xif) the number of connected components of IJj^]^ B{xi, 2u~^) with 
exactly I balls. Given rii, . . . , uk such that X]i=i ^"-^ ~ -^' ^^ define 

WK{ni, ...,nK)^ {XA- e T4^^, ^((xk) = n, for all I < K}. (20) 

Let us note that the subsets WK{ni, ■ ■ ■ ,nK), with J2i=i ^"-^ — ^' partition W^ . We then deal with two cases. 
1. If B{xi,2uj) n B{xj, 2uj) — for all i 7^ j < K i.e. x^ £ Wk{K, ■ ■ ■ ,0), the events considered in the 



right-hand side of (19 1 are independent 



2. If not, we are going to show that the contribution of such xa' in (14 1 is negligible 



More precisely, we write the integral (14) in the following way 



7^ / P (V^ < K,R{Cx^u{^^.}{x.)) < %) dxA = 7"^ / P (Vi < K,R{Cx,u{.^}{x^)) < u^) d^K 

+ 1" I P (Vi < K,R{Cx-,xM^^}{xi)) < u-,) d^K- (21) 

JW'^-WKiKfl,...fi) 



Step 2.1 (Case of disjoint balls) For all xa" = {xi,. . . ,xk) £ Wk{K,Q, ■ ■ ■ ,0), we obtain from (19 1 and (18 1 

K 

7^P (V^ < K,R{Cx,u{.^}{xi)) < u,) ^l[jP{R{Cx.,u{.M^r)) < u,) ^^ {t^+^)'' . (22) 

i— 1 



Moreover, \dK{WK{K, 0, . . . , 0)) — > 1. This shows that 



l"" / P (V* < K,R{Cx,^{^^}{xi)) < u^) d^K -^ (t'^+i)^. (23) 

JWK{K,0,...fi) ^^°° 

Step 2.2 (Case of non disjoint balls) In this step, we show that the second integral in the right-hand side of 



( 21 1 converges to 0. In particular, we study the limit behaviour of the integrand of ( 14 1 for all xa- = (xi , . . . , xk) £ 
WK^ni, . . . jUk) with (ni,...,nA) 7^ (-R', 0, . . . , 0). The number of points of X^ D [J^^^^ B{xi,2uy) is Poisson 
distributed of mean ^Xd ( Ui=i B{xi, 2u^) j . From ( [T9| , we deduce that 

7^P (Vi < K,R{Cx^u{^^}{x,)) < u.,) 

V LAl^ li^^-.A.(U.^,B(..2„-,)) xp^(^^^...^^^). (24) 



' ^ fc! 

A:=0 



The term Pk{xi, . . . , xk) denotes the probabiHty to cover the spheres S{xi^ u^), i = \ . . . K, with the spherical caps 
{ J24 ■{xi),i ^ j < K} and {^y„ (xi),m < fc} , defined in ( 16 1 , where j/i , . . . , y^ are k independent points which are 
uniformly distributed in IJj^]^ B(xi, 2u^). This probability satisfies the following property: 

Lemma 3. Let x/^ ~ (xi, . . . , xk) S Wj^ini, . . . , n^) and 



K 



iV = 5](d+l) 



ni. 



(25) 



1=1 
Then, for all /c < TV 



Pfc(xi,...,x_R-) = 0. (26) 

The proof of Lemma p] is postponed to the appendix. From (24), (26 1 and the trivial inequalities < 
Pk{xi, . . . ,xk) < 1 and Xd [[J-^i B{xi,2u~^)] < k2'^K,du'^, we deduce that there exists a constant c, depending 
on K, such that 

7^^P (V. < K, i?(Cx,u{x..}(-^)) < -7) < 7^ E ^^"'"V"'^' -00 ^ ■ ^''(^O'' 



0(7) 



where (/)(7) ~ ^"(7) means yP{ — > 1. Using (11), (25) and the fact that K = y2l=^ ^""-i^ ^^ obtain for 7 large 

7->oo vl7J -y— ^00 ' ' ' ' 



-K 



enough 

7^P (V. < K,R{Cx,u{.^}{x.d) <u,)<c •[] 7(^-1)"' . 



K 



(27) 



1=2 



K (ei+2)(i-l) ^ 



Moreover, using the fact that AdK(WV("^i, ■ ■ ■ , nx)) < c ■ ni^2("7)^'^^^"' = c ' 11^2 7^' '^+^ '"' and ([27]), we 



get 



7^ / P (Vz < if, i?(Cx,u{xK}(^0) < %) rfxK 

Jvy^-WK(-ftr,o,...,o) 



E 



„K 



K 



7- / P(Vz<if,i?(C;f^u{xK}(2^.))<%)rfxK<c-En7"^' 

WK[ni,...,nK) 1—2 



(28) 



The sum above runs over all the if -tuples (ni, . . . , tik) such that X]i=i ^'^i — ^ and rii ^ K. Since (ni, . . . , tik) 7^ 
(if, 0, . . . , 0), there exists I > 2 such that n; 7^ 0. Consequently, we get from (28 1 



,K 



1" I r (Vi < if, R{Cx,u{^^}{xi)) < u,) d^K - O (7-1/(^+1)) 

IW'<-WK{Kfi,...fi) ^ ^ 

4>h) 



(29) 



where ^(7) — 0{Tp{l)) means that ;?p^ is bounded 



Conclusion From ( 23 1 and ( 29 ) , we deduce that for all if > 1 

7^ / P (V* < if,i?(Cx,u{x^}(a;.)) < %) d^K -^ (t^+i)^. 
We then apply Lemma [2] with A= R, to conclude that 
P (i?mi„(7) > W7) -^ e-*'^\ 

7— >-oo 
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The cell which minimizes the circumscribed radius is asymptotically a simplex. To show it, we denote by 

where Fd~i{Cx .{x)) is the number of hyperfaces of Cx {x). The order of convergence of i?I„in(7) is greater than 
u^ according to the following proposition. 

Proposition 4. Let X-^ be a Poisson point process of intensity 7 and W a convex body of volume 1. Then, for all 
t> 0, 



P(a2«:.7^'+''/^'+''<i„(7)>i 



— > 1. 

7— >oo 



Proof of Proposition 4l We apply Lemma l2| to f{Cx-,{x)) — R{Cx^{x)) and A — [d + 2, 00). We then study the 
finite dimensional distributions i.e. 

7^ / P (Vi < K,R{Cx^u{^,M)) < "7'^d-i(Cx,u{x^}(a^.)) > d + 2) dx^ (30) 



for all iiT > 1. When X = 1, the integrand of (30 1 is 



7P(i?(Cx,u{0}(0)) < u^,Fd-,{Cx^u{o}m >d+2)< jP {R{Cx.,u{o}m < u^,#{X^ D B{0,2u^)) > d+2) 

•"^-^ k\ 7-5-00 

k=d+2 

We deduce that 7/^^ P {R{Cx.,u{x}{x)) < Uj,Fd^i{Cx,,u{x}{x)) > d+2) converges to 0. More generally, for all 
iiT > 1, we get 

7^^ / P (V* < K,R{Cx^u{^^M)) < u^^Fd-iiCx,u{.^}ix,)) >d + 2)d^K ^ 0. (31) 

JWK(K,o,...,a) '<^^ 



Moreover, from (29) 



7^^ / P(Vi < K,R{Cx,^{^^}{x,)) < ?^7'^d-i(^x,u{xK}(^«)) >d+2)d^K ^ 0. (32) 

From (31 ), (p2|) and Lemma [2] applied to A = [d+2, 00), we get 



' (i?:„i„(7) > U7) ^ 1. 

7—^00 



n 



Corollary 1. Let X^ be a Poisson point process of intensity 7 and W a convex body of volume 1. Then 
P (Vx e X^,R{Cx,{x)) =. Rn,inh) =^ i^d-i(Cx,(x)) = d+l) -^ 1. 

7— >oo 

Proposition El implies Corollary fTl but does not provide the exact order of R'nunil)- Nevertheless, when d = 2, 
it can be made explicit. The key idea is contained in Lemma Hand cannot unfortunately be extended to higher 
dimensions. 
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Proposition 5. Let Xj be a Poisson point process of intensity 7 and W a convex body of volume 1 in R^. Then, 
for all t > 0, 



where cJ^ ^^ defined in ( 35 1 



Proof of Proposition [5[ Let i > be fixed and let us denote by 



4(i) = [d^'.-'^c'i't) 



1/2 



(33) 



where a^ is specified in (35 1. As in the proof of (2d I, we interpret the distribution function of -Rmin(7) ^^ 
a covering probability of the circle. Let /i^ be the probability that S'(0,u' ) is covered with the circular caps 
{j2^^^(0),m, < fc} where j/i, . . . , y^ are k independent points which are uniformly distributed in i?(0,2u' ) and such 
that'>i(qo}u{y.}(0))>4i.e. 



00 ^ 

(i?(Cx,u{o}(0)) < <,Fi(Cx,u{o}(0)) > 4) = ^ -(4^7w;2)fcg-4.7<^^. 



(34) 



/c=4 



The constant Og is defined as 



32 



1/4 
M4 I > 0. 



(35) 



We are going to apply Lemma [21 to the event A = [4, cxj) replacing u^ by u' . To do it, we need to get the limit 
behaviour of 



7^ / P (Vz < X,i?(Cx^u{x^}(2;.)) < <,Fi(Cx,u{xk}(2^.)) > 4) d^K 



(36) 



for all K>1. 

When K ^ I, from ^ and (|3|), we deduce that 7 J,^ P {R{Cx^ij{x}{x)) < u'^, Fd^i{Cx^u{x}{x))) > 4) dx 
converges to t*. More generally, for all -ftT > 1, 



7^^ / P (V* < K,R{Cx.,u{^^}{x^)) < u;,Fi(Cx,u{xK}(a:,)) > 4) dx^ -^ t^^ 

JWK{Kfi,...,Vi) "^"^ 



7— >oo 



(37) 



Otherwise, for all x/f G Wxini, ■ ■ ■ , nx) with (ni, . . . , uk) 7^ {K, 0, . . . , 0), the integrand of (36 1 is 



7^P (Vi < K,R{Cx,^{^^}{xi)) < u'^,F,{Cx,u{^^}{x,)) > 4) 



7"E 



fc=0 



V LA1_ Ii^e-^'HU..^(^-M x^fe(xi,...,x^). (38) 



The term fik{xi, . . . ,xk) denotes the probability that S{xi,u') is covered with the spherical caps {s!/xj{xi),i 7^ 
j < K} and {.s!/y^(xi),m < k} where yi,...,yk are k independent points which are uniformly distributed in 
[J^^^B{xi,2u' ) and such that Fi{Cx u{xK}ixi)) > 4 for all i < K. This probability satisfies the following 
property: 
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Lemma 4. Let x/f = (xi, . . . ,xk) € Wxini, . . . ^uk) C R^ and 

K 

iV' = 4ni + 4^2 + ^ ini . (39) 

Then, for all k < N' 

lik{xi,...,XK) ^0- (40) 



The proof of Lemma [i] is postponed to the appendix. From ( 38 ) , ( 40 1 and ( 39 1 , we deduce for 7 large enough 
that 



K 

41-3 , 



7^P (V^ < i^, R{Cx,u{.^}ix,)) < u\, ^i(Cx,u{x^} (2;.)) > 4) < c • 7'^(70'^' = c • 7"^ H ^"^"' 



;=3 



Moreover, \2K{WK{ni, . . . .uk)) < c • n^2("7 )^'"^^"' = c • 7-3"^ n/=3 7^^"'- This shows that 



l"" / P (Vz < if,i?(Cx,u{xK}(a^.)) < n^,F^{Cx-,v^{^^}{x^)) > 4) dx^ - O (7^1/*) . (41) 

JH'^-H'k(-R',0,...,0) ^ -^ 

From ([37]) , ([4l|) and Lemma [2] we get 



i?Ln(7) > < ^ e 

' 7—^00 



-f 



n 



We conclude the section with a quick sketch of proof for ( 2a ) 



Proof of (2a I. We notice that 



''max(7) = max riCx {x)) — - max min d(x,y). 

The behaviour of the maximum of nearest neighbor distances was studied by Henze in Theorem 1 of (T3] when the 
input is a binomial process. His result did not include the contribution of boundary effects and is consequently 
limited to the set of points in VF 5(0, w^). With Lemma Inland proceeding along the same lines as in the proof 



of ( 2d I , we are able to show the convergence in distribution of the maximal inradius of Voronoi tessellation when 



the input is a Poisson point process in VF. D 



4 Proof of (2c), consequence on Poisson- Voronoi approximation 



Proof of (|2cl). First, we notice that 



-Rmax(7) = max R(CxAx))= max max d(x,y). 

xf^x^nw ^ x£X^nWyi^Cx^(x) 

In order to avoid boundary effects, we start by studying an intermediary radius -Rj„ax(7) defined as 

i?' ,^(7) = max max dix.y). 

xex-,,cx.,(x)nw^0yecx^{x)nw 

In a first step, we provide the asymptotic behaviour of R'^^^{'^)- Secondly, we study the effects of Voronoi cells 
astride W and W^. 
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step 1 The distribution function of i?max(7) can be interpreted as a covering probability. Indeed, if we denote 
by 

/ 1 1 \^^'' 

u^^u^it) = 1+ log(ai7(log7)''-i) (42) 

where 

and t is a fixed parameter, we have 

<iax(7) < % <^=^ Vx e X^, s.t. Cx.,(x) n P^ 7^ 0,Vy e Cx^{x) n M^,d(a;,2;) < u^ 
-<=> Vy e W, 3a; e X^, (i(x, y) < u^ 
<=> {5(2;,^^), a; e X-y} covers PF. 

We have to deal with the probability to cover a region with a large number of balls having a small radius when 
7 — > 00. Asymptotics of such covering probabilities have been studied by Janson. We apply Lemma 7.3 of |17j 
which is rewritten in our particular framework. Actually, Lemma 7.3 of 17J investigates covering with copies of a 
general convex body and requires conditions which are clearly satisfied in the case of the ball (see Lemmas 5.2, 5.4 
and (9.24) therein). 

Lemma 5. (Janson) Let W he a bounded subset of R'' such that Xd{dW) = and Xj a Poisson point process of 
intensity 7. Let R a random variable such that E[R] > and E[i?'^+'^] for some e > 0. We denote by a{B{0, R)) = 
aiE[R''-^]'^E[R'^]-'^'^-'^'> . If a = a{j) is a function such that 0(7) — > and 

then 

P{{B{x,R),x£X^} covers W^) — > e^"''. 

7— >oo 

Taking a = u^, R — I, Xd{W) — 1 and noting that E[Xd{aB{0,R))] — KdU^ and a{B{0,R)) = ai, we check 

easily (44 1. From Lemma ^ we deduce that P {{B{x,u^),x e X^} covers W) converges to e^'^ . Hence, for all 

te R, 

lim P(i?:„,,(7)<«^) = e-^^'. (45) 

7— J'OO 

Step 2 Taking f{Cx^{x)) = Kd{i-na.yiy^Cx-,ix)nw d{x,y)y\ a^ = 7, 5-, = 7log (ai7(log7)''"i) and Y a Gumbel 
distribution (i.e. P(F <t) — e"'^ ,f e R), one can check condition ^ with k — d. From ( [is] ) and Proposition^ 
we deduce that Eirnax-x^x^nw '^^s^'^yeCx (x)r\wd{x,y) < Uj) converges to e"*^ for all t G R. Using the fact that, 
on the event A.y (given in Lemma fll , 

max max d(x,y) < max max d(x.y) < max max d(x,y) 

xex-,nwvecx^(x)nw xex-,nw yeCx^ix) xex^nWi+i^ yeCx-,{x)nWi+i^ 

and proceeding along the same lines as in the proof of Proposition l3] we get 

P (i?max(7) < U7) -^ e-^'' . (46) 



7—^00 
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We can note that the asymptotic behaviour of i?max(7) gives an interpretation of Lemma 7.3 in |17j . Indeed, 



( 46 1 shows that the Gumbel distribution which appears as a hmit probabihty of a covering is actually the limit 
distribution of a maximum. 

We now apply this convergence result to the so-called Poisson-Voronoi approximation defined as 

It consists in discretizing a given convex window W with a finite union of convex polyhedra. This approximation 
has various applications such as image analysis (reconstructing an image from its intersection with a Poisson point 
process, see [19^ ) or quantization (see chapter 9 of [10 J. Estimates of the first two moments of the symmetric 
difi^erence between the convex body and its approximation are given in |14j and extended to higher moments in 
[33]. To the best of our knowledge, the convergence of Yx {W) to W in the sense of Hausdorff distance, denoted 
by dni', ■), has not been investigated. Corollary p] addresses that question with an assumption on the regularity of 
W which is in the spirit of the n- regularity (see Definition 3 in [7]). 

Corollary 2. Let us assume that there exists a > such that, for v small enough and for all y E W, 

XdiBiy,v)nW)>aXdiB{y,v)). (47) 

Then 

p(dH{W,yx,{W))<{c{a)j-'log{aij{log^)''-')y^'') -^ 1 (48) 

where c(a) = k^^ + 2'^K^^a~^. 

Proof of Corollary |2} Let us denote by 

v^ = (c(«)7-ilog(ai7(log7)'^-i))'/'. (49) 

First, we show that max^g ^^ (vi') d{y, W) < v^ with high probability. For all t € R, using the fact that m^ < v^ for 



7 large enough, where Uj = Uy(t) is given in (42 1, we get 



P max d{y, W) < V^ ] > P (i?,„ax(7) < V^) > P (i?n.ax(7) < M7) . 

\yerx^{W) 'J 

From (46) and Proposition k5I the last term converges to e~^ as 7 goes to infinity. Taking i — >■ c», we get 

lim P I max d(y,W) <v^] > lim e"""' == 1. (50) 

7-+00 \ye'rx-,{W) 'J t^oo 

In a second step, we are going to show that maxygvK d{y, X^ D W) < v-y with high probability via the use of a 
covering of W by balls as in the proof of (pc|. Now, the convex body W is covered by ^ = O {v^'^) deterministic 



balls Bi, . . . , B^ with center in W and radius equal to v^/2. From H7j, (49 1 and the fact that #(5^ n {X^ n W)) 



7/ 

is Poisson distributed with mean jXd{Bi n W), we get for 7 large enough 

■J/ 



('max%,X^nVF) > v\ < P ( Q{#(B, n (X^ n I^)) = 0} j < ^e-T""''^"-/^)' < arS"'(log7)"^''"'^^. 
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Using the fact that .yV = 0{v '^) i.e. ^ — 0(7(log7) ^) according to (49), the right-hand side is O ((log7) ''). 
Hence 



I max(i(y, X^ n W) < w- 



\yew 



1. 



(51) 



Since dH{W,'fx-,X^)) < max -^ maxj^g ;^_^ j^-iy)d{y,W),Taayiy^iY d{y,X^ nW)>, we deduce from (50 1 and (51) that 
P{dH{W,yx,{W))<v^) -^ 1. 



n 



In [13], Hevehng and Reitzner obtain that the volume of the symmetric difference between W and 1^ (W) is 
of the order of 7^^/'*. The result above makes sense and could provide the right order of the HausdorfF distance. 
Obviously, the constant c{a) — k'^ + 2''^^ a~^ is not optimal. From Lemma 111 it would have been possible to get 
an upper-bound of the order of 7^(1^'^)/'^ but it is less precise than Corollary [2| 



5 Proof of (2b) 



Proof of (2b I. Let t > be fixed. We denote by Uj the following function: 

l/d 



,(t) ^ (2-('^-i).,-V-t) 



(52) 



We start by finding a different expression of r,„in(7) which does not rely on the Voronoi structure. Indeed, for all 
X G Xy n H^ we have 



r(Cy (.t)) = max{r > 0,B{x,r) C Cy (.x)} = - min d{x,y). 

Hence, rmin(7) can be rewritten as 

rmin(7) = o ™^ d{x,y). 

2 (x,y)^e(X-,nw)-KX., 



(53) 



The equality (53 1 implies that the problem is reduced to a study of inter-point distance. Such study is well known 



for a binomial process X*^"' of intensity n in W . In particular, Jammalamadaka and Janson (see [IS], §4) have 
shown that for all t > 0, 



P ('^min.n > ^«) 



where r'^^-^ „ is defined as 



1 



(54) 



2 (x,y)^eX<")xX(") 



d{x,y) 



and Un given in (52). In a first elementary step, we extend the limit to a Poisson point process. Our main 



contribution is then to compare the obtained limit with ri„in(7) by dealing with boundary effects. In particular, 
our study provides a far more accurate estimate of the contribution of boundary cells (see ( 66 1) than what we could 
have deduced from Proposition [3] 
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step 1 We extend (54 1 to a Poisson point process. We define 

rLnil)^:^, , min d{x,y). 

Let us note that for all < a < /3 < 1, and for all n g {0, 1, 2, . . .}, |n — 7I < 7" 
enough. Consequently, since u-^ is non-increasing in 7, we have for 7 large enough 



(55) 
1^ ~ 7I ^ n^ for 7 large 



|P(C„(7) > U7) - e"*| < E IPKnin.n > %) " ^^l ^(#(^7 01^) = «) 

< J2 max {|P(r:„i„^„ > u„_„,) - e-*| , |P(r;„i„,„ > u„+„,) - e"*]} P(#(X^ r]W)^n) 

\n — |'I<7° 

+ P(|#(X^nW^)-7l >7"). (56) 



The second term of (56 1 converges to thanks to a concentration inequality for Poisson variables (see e.g. Lemma 1.4 
in [SI])- The first term is lower than max„>..y_^Q max 1 1 P(rj'j^;j^ „ > w„_„,8 (i)) — e^* | , | P{r[-^^^^^ „ > u„^„fi (t)) — e~* | } 
which tends to according to (54 1. This shows that, for all t > 0, 



lim P(Ci„(7) >M^) = e *. 



7— foo 



(57) 



Step 2 We show that r„^ij^{'^) — ^^^^^^(7) with probability of order of 0(7^') with e G (0, |). Indeed, the random 
variables rniin(7) and r'^-ninil), defined in (ISSl) and (55 1, are equal if and only if no point of X^ n W^ falls into the 
union of the balls B{x, 2r^i,j(7)) for x & X^nW such that d{x, dW) < 2r^ji„(7) i.e. 



/ / 



P(nni„(7) ^ r'„^nAl)) = P 



# 



V 



\ \ 



X^nW^n U B(x,2r^„(7)) 



V 



xex^nw, 

d(x,aW)<2r' . (t) 



^0 



/ 



/ 



< E 



J2 #(X^nW^^ni?(x,2rLn(7))) 

ix^n' 

f)<2r 

From Shvnyak-Mecke formula (see e.g. Corollary 3.2.3 of [27]), we get 



xex^nw, 

d{j:,aW)<2r' . (-,) 



(58) 



53 #(X^nM/^nS(x,2r^„(7))) 



xex^nw, 



7E 



where r^,Zil) = ^t^^t^(x' ,y)^e{X-,u{x}nw)'^ d{x',y) for all x G X-^DW. Noting that r'^J^j) < r'^inil), we then 



# (x, n w^^ n i?(x, 24^(7))) i,(,,,^)<2.:^^.(,) 



da; 



obtain 



53 #{X^nW^nB{x,2r[^,^{^))) 



d{x,0W)<2r' . {-)■) 



< 7E [# {X^ nW^n B{x, 2r'^,^{j))) td(x,dw)<2rUJi)] dx 



w 



w 



7E [E [# {X^ C^W'C^ B{x, 2r[^,^{-f))) \X^ n W] ldix,aw)<2r:^,J^)] dx. (59) 
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Since # {X^ n VF^^ n B{x,2r[^i^^{'^))) is Poisson distributed, we get 

7E [E [# {X^ nW^n B{x, 2r^„(7))) 1^7 ^ W] ld(.,8M/)<2r;„,„(7)] 

= 7^E [Ad(M^''nS(a;,2r^i„(7)))lrf(^,aH')<2r;^._^(7)] < 2'^K<i •7^E Cin(7) lrf(i:,9H')<2r;,i„(7) ■ (60) 



Using (58 1, (59), (|60|) and Fubini's theorem, we obtain 

Knin(7) / ^d{x.,dW)<2r'^^J-()dx 



P(n„in(7) ^ '^mi„(7)) < Z^Atd • 7'E 



IV 



<c-j'E[C+\j)]. (61) 



The last inequahty comes from Steiner formula (see (14.5) in [37]) and c denotes a constant depending on W. 
Hence, to show that 

P (rnunil) ^ r'^inil)) ^ (62) 

7— >oo 



we have to find some upper-bound of 7 E[r^j^ (7)]. We know, from (57l and (52 1 , that 7 ?'i„i„(7) tends to 
in distribution but it does not imply ( 62 1 . Lemma p^ below provides an estimate of the deviation probabilities of 
7^ ''mill (7) when the window W is a cube. 

Lemma 6. Let C be a cube of side M and Xj a Poisson point process of intensity 7. Let us denote by 

Ci„|c(7)= ^ . , ™™n^,.^(^'2^)- 
Then, for all u < min{|;Afd^/^, id^/^7"-^/''}, there exists a constant c{M) such that 
P('';i„|c(7) >u)< e-(*'^)-="^ 



Proof of Lemma [6} Let u < minjiAfdi/^, ^^1/2^-1/^} be fixed. 

We subdivide the cube C = [0, M]'^ into a set of ,yV subcubes Ci, . . . , C^ of equal size c with c = 2d~^^'^u and 
^ = [Mc~^J . Since diam(Ci) = 2w for each i < ^, we obtain 



■-^ ^ .^ 



nr'^in\ch) >u)<P mmC^nX,) < l}j = (e-^^''(l + 7c'^))^ 
Replacing c'' by 2''-d~'^/'^u'^ and J/' by [2~^Md^/^w~^J'* we obtain the following inequality: 

P(Cin|c(7) > ") < eL2-^^^''^^^"-^j'('°s(l+72^^-/^«^)-72^.-/^n^) 

Since -f2'^d-'^/^u'^ < 1 and 2M-^d-^/^u < 5, we have log(l + j2'^d-'^^'^u'^) - -f2'^d-'^/^u'^ < -\2'^'^d''^-i^u^'^ and 
[2-iA/di/2u-ij'^ > (2-iAfdi/2^-i _ ;L)d > 2-2dAf'*d'^/2^-'*. Hence 

P(Ci„|c(7) > ") < e-^'^-^'^^^^V"'' ^ e-(*^)-^"^ 
where c(M) = id^''/^^/''. D 
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Now, we can derive an upper-bound of 7^E[r^j^ (7)]. Indeed, since W has non-empty interior, there exists a 
cube C of side M included in W. Using the fact that #(X-, n C) > 2 ==> r[^-^ia{l) < '"minic^'^)' ^® S^* 

/'diain(iy) 

7'EKt„'(7)] - 7' / P(r-:^[i;\7) > s)ds 

Jo 

< dmmiWhMMX, nC)<l)+j' nCtU^) > s)ds. (63) 

Jo 



The first term of the right-hand side of (63) is decreasing exponentially fast to since #(X^ n C) is Poisson 



distributed of mean 7M'*. For the second term, let us consider a fixed e in (0, |). Then 



7- -f MdV2^2p (,;^^^|^(^) > ^-(2-H.)/(.+i)) . 



I ■ I \ mill I o V ' ^ — ' / 

(64) 

Since e > 0, we have 7~(2+€)/(d-i-i) <- minjiMd^/^, i(i^/^7~^/''} for 7 large enough. Hence, from Lemma 6 applied 
to u := 7-(2+<^)/(<^+i)^ we deduce that for 7 large enough, 

7' / P(vtn[c(7) > -)d. < 7"' + MdV2^2^-c(M).<--)/(-^) ^ (g5) 

Jo 



Combining that argument with mlh, (^63k and (651, we deduce that 



The last term of the right-hand side of ( 65 ) converges exponentially fast to as 7 goes to infinity since e < 3 



P (rnUl) ^ r'^^inil)) ^ O (j-^) . (66) 



We then deduce from ( 57 1 and (|66| that 



|P (r^„(7) > u^) - e-*| < |P (r^j7) > "7) " e"*| + 2P (rminil) + r[^Al)) ^ 0. 

I 'II 'I 7— >-oo 
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Remark . The rate for the convergence in distribution of rmin(7) to the WeibuU distribution can be estimated. 
For instance, we can show that Theorem 2.1 in [38 implies the rate of convergence of 7'j„j,j(7). Another way to get 
it is to use Theorem 1 in 1 . We then obtain that there exists positive constants c{W) and r(VF) such that, for all 
e < |, t > and 7 > r(Ty)', 

|P (2'^-iK<i7'''mi„(7)' > - e"*| < c(W^)7"™"^'"^ 

The study of more extremes for general tessellations and their rates of convergence will be developed in a future 
paper. 

Appendix 

Proof of Lemma Isl Actually, we show the following deterministic result: let JC > 2, A: < A^, (xi, . . . , xk) & 
WK{ni,...,nK) with (ni, . . . ,n_R-) ^ {K,0,...,0) and {yi,...,yk) G \J-^^B{xi,2uj) such that {xKJUJyfc} are in 
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general position i.e. each subset of size n < d+ 1 is affinely independent (see [10]). Then there exists i < K such 
that sphere S{xi,u-y) is not covered by the induced spherical caps {s/^ {xi),i ^ j < K} U {s2/y^{xi),m < k}. 



K 



Indeed, from (25 1 there exists a connected component of IJj^]^ B{xi, 2u^) of size 1 <l < K, say IJj^^ -^(^i; 2u^) 



without loss of generality, such that Ni < d + 1 with 

Nl = ^{{Yk]C^[jB{x,,2uM (67) 

Since {x;} U {yAf, } arc in general position, the family {x;} is not included in the convex hull of {yat, }. In 
particular, there exists i < I such that Xi is not in the convex hull of {x/} U {yni} — {xi}- Since a Voronoi cell 
induced by a finite number of points is not bounded if and only if its nucleus is an extremal point of the polytope 
induced by the points, it implies that the circumscribed radius of C{x]}u{yjv }i^i) ^^ '^ot finite i.e. S{xi,u-y) is not 
covered. D 



Proof of Lemmal4| We show the following deterministic result: let iiT > 2, fc < N', {xi, . . . , xk) & WK{ni, ■ ■ ■ , tik) 
with {ni,...,npc) 7^ {K,0,...,0) and (j/i, . . . , 2/fc) € [j^^^B{xi,2uj) such that {x^} U {yj;} are in general posi- 
tion. Then there exists i < K such that either the sphere S{xi,u-y) is not covered by the induced spherical caps 
Wx,{xt),i^ j < K}(j{£/y^{xi),m < k} or Fi{C{^^}u{yk}ixi)) < 3- 



K 



Indeed, from (39 1 there exists a connected component of lJi=i B{xi, 2u-y) of size 1 <l < K, say (Jj^^ B{xi,2u^) 



without loss of generality, such that Ni < 4 if Z = 1 , 2 and iV; < 3 if Z > 3 where Ni is given in ( 67 1 . 

• If / = 1, either S{xi,u^) is covered or Fi(C{x^}u{yfc}(a;i)) = Fi{C{^^}yj{yj^^}{xi)) < 3 since Ni < 3. 

• If / > 3, from Lemma ^ there exists i < I such that S{xi,u^) is not covered. 

• If / = 2, we can assume that A'2 = 3. We have to prove that if y^ = {j/i, ?;2, Us} is a set of three points in 
B(xi,2u^) U B(x2,2u^), then the following properties 1 and 2 below cannot hold simultaneously. 

1. The circles S{xi, u^) and 5(x2, u^) are covered by the induced circular caps 

WxAx2),J^X2{.Xi)} U {^y„(xi),rn < 3}. 

2. The number of edges of the Voronoi cells satisfy Fi(C{a,j 3,, y3}(a;i)) > 4 and Fi{C{rj.^^^y^^{x2)) > 4. 

Let us assume that Properties 1 and 2 hold simultaneously. Let us denote by G the Delaunay graph associated 
to {a;i,X2, 2/1, 2/2,2/3}- Then G is a connected planar graph with 1; = 5 vertices and e edges. From Euler's 
formula on planar graphs, e < 3v — 6 i.e. 

e < 9. (68) 

From Property 1 and according to the proof of Lemma ^ a;i,a;2 are in the convex hull of {2/1,2/272/3} i-^- 
{2/1,2/2}, {2/1,2/3} and {2/2,2/3} are edges of the associated Delaunay triangulation. From Property 2, a;i,X2 
are connected to every point i.e. {xi,a;2}, {xi,yi}, {^1,2/2}, {xi^ys}, {x2,yi}, {3:^2,2/2} and {a;2,2/3} are also 



edges of the Delaunay triangulation. The total number of these edges is e = 10. This contradicts (68 1 



D 
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